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I.  The  current  density  distribution,  flow,  pressure,  and 
density  fields  are  calculated  for  a  (quasi- )__  steady  nozzle- 
type  MPD  thruster,  based  on  a  partially  two/dimensional  ap¬ 
proach,  assuming  a  fully  singly  ionized  plasma  with  an  iso¬ 
thermal  electron  and  an  adiabatic  ion  behavior.  In  addition, 
for  a  giv^n  current  rise,  the  current  contour  lines  and  the 
flow,  pressure,  and  density  fields  within  the  cylindrical 
flow  discharge  region  are  calculated  as  a  function  of  time. 
As  initial  conditions,  one  assumes  a  homogeneous  radial  cur¬ 
rent  within  a  confined  axial  electrode  region  and  a  super¬ 
sonic  homogeneous  axial  flow. 

o  II.  In  order  to  investigate  the  arc  stability,  a  unique  the¬ 
oretical  approach  is  presented  which  allows  determination  of 
the  stability  behavior  of  an  arbitrarily  shaped  current- 
carrying  plasma  channel  under  various  conditions.  Based  on 
this  novel  theory,  the  onset  phenomenon  observed  in  any  noz¬ 
zle-type  MPD  arc  jets  can  be  predicted;  calculated  and  meas¬ 
ured  results  agree  favorably.  Moreover,  this  stability 
theory  can  also  explain  the  very  strange  behavior  of  arc 
spots  and  the,  until  now,  seemingly  obscure  retrograde  and 
Robson  drift  motion  of  arcs  in  a  transverse  magnetic  field. 


Ill,  The  cathode  attachment  of  electric  arcs  are  theoreti¬ 
cally  and  experimentally  investigated.  The  erosion  effects 
and  the  surface  damage  due  to  spots  are  discussed;  they  stem 
from  evaporation,  droplet  ejection,  and/or  splashing.  The 
most  severe  damage  to  the  cathode  surface  occurs  when  many 
samll  spots  cluster  together  and  build  one  or  even  several 
larger  spots  which  now  cause  melting  and  splashing.  Erosion 
rate  measurement  on  different  surfaces  and  cathode  materials 
are  presented.  / 
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Introduction  and  Obiectives 


The  on-board  power  requirement  for  future  space  missions 
will  soon  grow  into  a  range  high  enough  such  that  magneto- 
plasmadynamic  propulsion  systems  become  attractive  for  many 
near-earth  applications,  for  which  modest  specific  impulse 
values  of  10  to  30  km/s  are  optimal.  However,  the  develop¬ 
ment  stages  of  the  different  possible  MPD  thrusters  and  sys¬ 
tems  are  insufficient  to  reliably  predict  performance,  life, 
weight  and  cost  potentials.  Moreover,  the  fundamental  under¬ 
standing  of  MPD  discharges  is  still  incomplete.  Therefore,  a 
detailed  investigation  of  the  arc  discharge  and  acceleration 
processes  within  a)  the  electrode  gap  and  b)  the  arc  elec¬ 
trode  interphase  region  of  those  devices  is  mandatory. 

The  research  work  covered  in  this  report  deals  with  both 
problem  areas  and  intends  at  first  to  give  proper  explana¬ 
tions  of  the  still  not  fully  understood  limiting  processes 
occuring  in  MPD  arc  propulsion  devices  and  eventually  to  de¬ 
velop  design  criteria  for  annular  self-field  MPD  arc 
thrusters . 

Objectives 

It  is  necessary  to  gain  an  understanding  of  the  fundamental 
processes  occuring  in  coaxial  arc  jet  and  MPD  arc  thrusters, 
in  order  to  obtain  reliable  design  criteria  and  to  predict 
thruster  performances. 

To  achieve  these  goals,  one  must 

-  investigate  the  arc  cathode  and  anode  attachments  under 
different  operating  modes  (experimental  and  analytical), 

-  develop  computational  methods  to  determine  the  flow  arc 
discharge  region  and  the  acceleration  mechanisms,  and 

-  analyze  the  arc  stability  under  various  operating  condi¬ 
tions  and  geometries  (onset) 

Based  on  these  objectives,  the  report  is  divided  into  three 
sections:  »flow  discharge  calculations, 

•stability  theory  and 
•electrode  investigations. 


2 .  Flow  Discharge  Calculations 

2 . 1  Nozzle-Type  MPD  Thruster  Investigations 
(Steady-State ) 

In  order  to  predict  the  overall  performance  of  a  nozzle-type 
self-field  MPD  thruster  (Figs.  1  and  2)  a  semi-two-d imen- 
sional  model  calculation  and  computer  program  has  been  de¬ 
veloped.  In  this  approach,  the  current  contour  lines  are 
calculated  in  r-  und  z-  dependence  while  the  flow  field  is 
taken  as  a  one-dimensional,  frictionless  expansion  flow  as¬ 
suming  isothermal  behavior  for  the  electron  component  and 
isentropic  expansion  for  the  heavy  particles  like  ions  and 
neutrals.  This  assumption  is  justified,  since  primarily  the 
electron  component  is  heated  by  Ohmic  energy  input,  and 
since  within  the  nozzle  expansion  flow  the  energy  exchange 
between  the  lightweight  electrons  and  the  heavy  particles 
like  ions  and  neutrals  is  fairly  weak.  Likewise,  reaction 
collisions  (ionization  and  recombination  processes,  etc.) 
are  neglected.  Hence,  the  model  approach  distinguishes  be¬ 
tween  electron  and  heavy  particle  temperatures  and  assumes  a 
frozen  flow  from  the  nozzle  throat  on  downstream. 

Within  the  nozzle  throat,  one  assumes  a  "fully  developed" 
arc  flow  regime,  i.e.  the  gas  and  plasma  velocity  has  only 
an  axial  component  and  the  arc  column  is  defined  by  only 
radial  heat  transfer  (and  possibly  by  radial  radiation) 
losses.  Moreover,  from  the  nozzle  throat  on  upstream,  one 
no  longer  distinguishes  between  electron  and  heavy  particle 
temperatures,  since  within  this  higher  pressure  regime  a 
tight  coupling  between  all  components  is  more  likely  to  be 
established  than  in  the  low  pressure  downstream  nozzle  re¬ 
gime.  In  spite  of  the  neglect  of  the  radial  dependency  of 
the  flow  and  species  temperatures,  however,  this  model  ap¬ 
proach  turns  out  to  be  an  effective  tool  for  predicting  the 
current  contour  lines  within  and  the  thrust  and  specific  im¬ 


pulse  of  a  nozzle-type  self-field  MPD  thruster. 


The  thrust  of  a  self-field  MPD  accelerator  is  composed  of 
two  effects;  first  the  sum  of  all  gas-dynamic  surface  forces 
on  the  inner  plenum  chamber  and  on  the  inner  nozzle  wall, 
and  second  the  sum  of  all  volume  effects  given  by  the  mag¬ 
netic  field  forces  which  act  between  the  current-carrying 
plasma  and  the  entire  power  circuit,  including  the  coaxial 
electrodes.  Hence  it  is 

Th  =  /  {pvv  +  P}*dA  -  f  j  x  §  dV  (1) 

As  V=Vi+V2 

=  c  ,p0A  /  j  x  B  dV  (la) 

V=Vi+V2 

where  (pvv  +  p}  is  the  gas  dynamic  force  per  unit  area  and 
As  represents  the  surface  of  all  internal  walls  like 
plenum  chamber  with  the  gas  inlet  cross  sections,  the  cath¬ 
ode  surface  and  the  inner  solid  wall  of  the  nozzle.  The 
second  term  is  the  sum  over  all  j  x  &  forces  which  act  be¬ 
tween  the  current-carrying  solid  leads  of  the  electric  power 
circuit  and  the  current-carrying  plasma;  in  the  second 
equation  (la)  po  is  the  chamber  pressure,  At  the  nozzle 
throat  area  and  cph  the  thrust  coefficient  which  matches 
the  gas  dynamic  effects.  It  is  by  definition 

c'rh  =  /  {pvv  +  P}*dX  (2) 


The  above  thrust  formula  eq.  (1)  can  be  rearranged  by  means 
of  the  equation  of  motion 

V  { pvv  +  P  }  =  3  x  B  (  3 ) 

which  is  integrated  over  the  volume  Vi  of  the  plenum  and 
nozzle  (see  Fig.  3).  With  the  theorem  of  Gauss  follows 


Powersupply 


Replacing  A/  {pvv  +  P}‘d^  -  ]  x  I  dV  in  the  formula  for 

s  l 

the  thrust  (1)  according  to  this  relation  (4)  one  obtains 


Ih  =  - 


{ pvv  +  P } *dX 


-  /  3  x  £  dV 

V2 


(5) 


or  in  our  one-dimensional  flow  approach,  by  neglecting  vis¬ 
cous  effects  and  by  taking  the  ambient  pressure  to  be  zero, 
the  absolute  thrust  is 


Th 
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I  /  ]  x  B  dV  j 
V2 


(6) 


where  now  m  is  the  mass  flow  rate,  vex  the  exit  velocity  and 
pex  the  average  pressure  at  the  nozzle  exit,  and  V2  is  the 
volume  of  the  current-carrying  "dome"  region  outside  the 
nozzle.  Hence,  instead  of  knowing  the  flow  and  pressure  on 
the  inner  solid  walls  of  the  thruster,  one  needs  only  to  de¬ 
termine  the  flow  and  pressure  on  the  exit  area  Ag  of  the 
nozzle  and  the  3  x  §  force  distribution  within  the  dome  re¬ 
gion.  Nevertheless,  the  3  x  §  distribution  should  be  known 
in  the  entire  current-carrying  plasma  region  of  volume  Vi 
and  V 2 ,  since  the  expansion  processes  within  the  nozzle  ( V 1 ) 
depend  on  the  axial  effects  of  the  3  x  S  forces.  This  is 
shown  in  Appendix  I.  One  obtains  for  the  thrust 
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and  hence  the  specific  impulse  by 


h 


where  v^  is  the  average  velocity  within  the  nozzle  throat, 
At  is  the  throat  area,  a  is  the  degree  of  ionization  and 
<  =  Cp/cv  the  quotient  of  the  specific  heats. 

The  calculation  of  the  thrust  and/or  the  specific  impulse 
and  the  current  contour  lines  follows  now  according  to  a 
stepwise  approximation  based  on  the  following  equations  and 
assumptions . 

Assumptions: 

•  steady  state  conditions 

•  rotational  symmetry 

•  singly  ionized,  quasi-neutral  plasma 

•  thermal  equilibrium  from  cathode  to  nozzle  throat--f rozen 
flow  from  there  on  downstream 

•  ambipolar  expansion  flow  without  friction 

(electrons:  isothermal;  ions:  adiabatic) 

•  within  nozzle  throat  fully  developed  flow  conditions  with 
equal  electron-  and  ion  temperatures 

•  electric  current  density  normal  to  electrode  surface 

Basic  Equations: 

Maxwell's  eqs.:  V  x  B  =  nol  a) 

V  x  E  =  -  =  0  b) 

>  n  (9) 
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Geometry  of  Thruster 


A(  z) 


Input  data: 


•  Mass  Flow  Rate  m 

•  Electric  Current  I 


Result : 


•  Current  Contour  Lines 

•  Thrust,  Th 

•  Specific  Impulse,  ISp 

The  stepwise  solution  proceeds  now  in  the  following  manner. 

First,  one  rewrites  Ohm's  Law  (10)  by  means  of  Maxwell's 
equations  (9a-d)  and  obtains  a  vector  equation  for  the  mag¬ 
netic  induction  vector  §  in  the  form 

-  =  —  curl  (—  curl  b)  -  curl  [v  x  B] 

dt  no  '•a  J 

+  —  curl  [3  (curl  B)  x  B]  =  0  (21) 

Ho 

Based  on  the  rotational  symmetry,  the  zero  azimuthal  current 
density  and  the  one-dimensional  flow  in  z-direction  along 
the  nozzle  axis,  one  obtains 
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and  for  the  quantity 
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follows  from  eq.  (19)  the  elliptic,  partial  differential 
equation  of  2nd  order  in  the  form 


d2¥  +  d2¥ 
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4 

dz 


(24) 


with  3 


CJT 

B  a 


(25) 


The  function  y(r,z)  -  const,  represents  now  a  current  con¬ 
tour  line,  since  B  =  B@  ~  I(r)/r  and  I(r)  is  the  electric 
current  carried  through  a  cross  sectional  area  of  nr2.  The 

proper  boundary  conditions  for  ¥  follow  from  the  geometry  of 

the  thruster  walls  and  electrodes. 

Step  1 :  By  setting  the  integral 

I  (7-K1  X  B)  A  dz 
2t 

in  the  expansion  relation  eq.  (18)  to  zero,  one  ob¬ 
tains  together  with  eqs .  (19)  and  (20)  a  first  ap¬ 
proximation  of  v  =  v(z),  p  =  p(z),  p  =  p(z)  and 

ni  =  ne  =  ne(z)  within  the  nozzle.  Outside 
the  nozzle  one  assumes  two  possible  flow  condi- 
t ions : 

Case  a:  an  extended  conical  plume  (vr  =  0)  and 
Case  b:  a  hyperbolic  expansion  flow  (vr  *  0) 


-li¬ 
as  shown  in  Fig.  4  and  Fig.  5,  respectively.  In 
both  cases  the  expansion  process  is  taken  adiabati- 
cally  for  electrons  and  ions  outside  the  nozzle. 

With  the  knowledge  of  the  flow  and  plasma  proper¬ 
ties  (in  first  approximation),  a,  3  and  v  (or  vz, 
vr  for  case  b  outside  the  nozzle)  are  also  known, 
as  well  as  their  derivations.  With  these  quantities 
one  now  solves  the  differential  equation  (24)  by 
means  of  a  finite  differential  method  (Gauss- 
Seidel)  for  the  current  contour  lines  in  first  ap¬ 
proximation.  In  Fig.  6abcd  and  Fig.  7abcd,  the  cur¬ 
rent  contour  lines  are  shown  for  Case  a  and  Case  b, 
respectively,  assuming  a  fully  ionized  argon 
plasma.  According  to  eq.  (9a),  one  also  obtains 
from  Y  the  current  density  and  the  j  x  §  force  con¬ 
figuration  within  the  discharge. 

Step  2:  One  solves  the  expansion  relation  eq.  (18)  numeri¬ 
cally  by  taking  ( j  x  §) z  within  the  integral  from 
Step  1.  With  that  the  velocity  v(z)  and  the  quanti¬ 
ties  p(z),  p(z),  ne(z)  are  calculated  in  second 

order  approximation  and  from  that  the  new  quanti¬ 
ties  for  a,  8  and  v  with  their  derivations. 

With  these  new  quantities,  one  now  solves  the  dif¬ 
ferential  equation  (24)  for  the  current  contour 
lines  by  means  of  the  same  numerical  method  as  in 
Step  1.  In  Fig.  8abcd,  the  current  contour  lines 
are  shown  for  four  current  settings  by  taking  a 
conical  plume  (Case  a)  outside  the  nozzle. 

Step  3:  By  taking  the  new  (  J  x  §)  configuration  from  Step 
2,  one  proceeds  now  in  the  same  manner  as  shown  in 
the  step  before. 

In  Fig.  9,  the  current  contour  lines  are  again 
shown  for  a  current  of  4000  A  and  compared  with  the 
results  of  Steps  1  and  2.  Based  on  these  results. 


4  jet  plume  conditions:  conical  expansion  flow 

with  constant  density  lines. 
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Fig  6d  Current  Contour  Lines 
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one  may  conclude  that  after  two  steps  the  current 
contour  lines  do  not  change  anymore  and  even  a 
first  step  approach  gives  good  results.  This  fact 
is  the  more  valid  since  for  lower  current  rates  the 
deviation  is  even  smaller. 

Compared  to  this  minor  deviation  between  the  results  of 
Steps  1  and  2,  the  deviations  between  the  contour  lines  of 
case  (a)  with  a  conical  plume  and  case  (b)  with  the  hyper¬ 
bolic  expansion  (see  also  Figs.  6abc  and  7abc)  outside  the 
nozzle  are  much  more  pronounced.  Nevertheless,  for  the  low 
current  of  1  300  A  both  contour  line  patterns  agree  quite 
well  with  each  other.  In  Figs.  10a  and  10b,  this  calculated 
current  pattern  is  compared  with  the  experimentally  measured 
contour  lines  of  the  same  current  and  show  good  agreement 
for  the  10%  and  even  20%  current  contour  lines.  The  fact 
that  in  the  measured  case  (Fig.  10b)  the  anode  attachment  is 
more  concentrated  downstream  to  the  nozzle  end  than  in  the 
calculated  case  (Fig.  10a)  may  be  attributed  to  the  theore¬ 
tical  model  which  neglects  the  boundary  layer  along  the  noz¬ 
zle  and  anode  walls. 

One  is  now  able  to  calculate  the  thrust  and  the  specific  im¬ 
pulse  according  to  the  known  j  x  ^  force  configuration  and 
the  flow  conditions  after  eq.  (7)  and  eq.  (8),  respec¬ 
tively.  The  results  of  these  calculations  are  shown  by  the 
upper  curves  in  Figs.  11  and  12.  Compared  with  the  measured 
values  they  are  by  a  factor  of  about  two  larger  for  small 
currents  (1000-2000  A)  and  approach  the  experimental  values 
within  less  than  20%  at  a  current  above  4000  A. 


By  properly  matching  the  thrust  coefficient  c^h  and  assum¬ 
ing  an  isothermal  behavior  for  the  electrons  and  heavy  par¬ 
ticles  from  the  chamber  on  downstream  to  the  nozzle  throat, 
one  may  calculate  the  thrust  according  to  eq .  (la).  The 
chamber  pressure  pa  can  now  be  related  to  the  nozzle  throat 
conditions  by 


a)  Calculated  and  b)  measured  current  contour  lines 
of  an  MPD  discharge  in  a  nozzle-type  plaana  thruster. 


Fig.  11  Comparison  between  calculated  and  measured  thrust 
of  the  Stuttgart  self-field  MPD  thruster  with  argon  as  pro¬ 
pellant. 
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(see  Appendix  3).  If  one  takes  c-ph  =  1,  the  thrust  follows 
by 

Th  =  1.65  mv  +  I  /  ]  x  §  dV  I  (27) 

V=Vi+V2 

The  result  based  on  this  formula  is  also  shown  as  the  middle 
line  in  Figs.  11  and  12. 

As  expected,  they  show  excellent  agreement  with  the  experi¬ 
mental  measurements  in  the  upper  current  range  (I  >  3000  A), 
and  even  down  to  1300  A  the  deviation  is  less  than  30%.  The 
pure  electromagnetic  thrust  and  the  corresponding  specific 
impulse  curve  are  given  by  the  lower  curves  in  Figs.  11  and 
12.  One  can  see  that  the  electromagnetic  thrust  reaches 
about  50%  of  the  total  thrust  at  a  current  of  about  4500  A. 

In  order  to  investigate  the  radial  effects  of  the  3  x  3 
forces  on  the  density  distribution  of  electrons  and  ions, 
one  assumes  a  balance  between  the  radial  [3  x  £ ] ^  and  pres¬ 
sure  forces.  The  results  of  these  calculations  are  shown  in 
Fig.  13ab  and  Fig.  14ab  for  the  two  cases  (a)  and  (b),  re¬ 
spectively.  In  these  model  calculations  the  expansion  proc¬ 
ess  is  taken  similarly  as  in  the  examples  before;  however, 
now  it  is  based  on  an  average  temperature  and  density  across 
a  nozzle  cross  section  or  across  the  ellipsoid  area  in  case 
(b)  for  the  flow  outside  the  nozzle. 
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(Time  Dependent  Calculation) 


In  addition  to  the  nozzle-type  MPD  arc,  a  detailed  numerical 
code  has  been  developed  which  allows  the  calculation  of  the 
plasma  properties  within  a  cylindrical  MPD  arc  as  shown  in 
Fig.  15.  This  experimental  thruster  consists  of  a  thoriated 
tungsten  rod  in  the  center  and  water-cooled  ring  segments, 
the  last  three  of  which  can  be  either  individually  or  all 
three  together  used  as  the  anode.  The  propellant  argon  (or  a 
hydrogen/nitrogen  mixture)  is  fed  through  a  small  opening  in 
the  backplate  and  enters  somewhat  downstream  of  the  flow 
discharge  region.  There  it  is  heated  up  and  axially  accel¬ 
erated  due  to  the  self-induced  magnetic  field  effects.  The 
thruster  can  be  run  in  pulsed  mode  or  continuously  in  steady 
state  mode. 

The  geometrical  configuration  on  which  the  computation  is 
based  is  shown  in  Fig.  16.  The  radii  of  cathode  and  anode 
are  taken  according  to  the  dimensions  of  the  real  device; 
however,  the  length  of  the  cathode  is  assumed  to  be  of  the 
same  length  as  the  anode  and  matches  the  length  of  the  three 
last  segments  (Anode  1,  2  and  3).  The  calculations  are  taken 
time-dependent  and  proceed  in  the  following  way: 

The  initial  conditions  are  the  total  mass  flow  mo,  the  tem¬ 
perature  To »  and  the  Mach  number  Mao;  likewise  the  current 
distribution  at  the  starting  time,  t  =  0  is  assumed  radially 
parallel  and  homogeneous  from  the  anode  to  the  cathode  (see 
dotted  line  in  Fig.  16).  The  current  rise  is  taken  according 
to  a  sine  square  function  (~  s  in  2  [  n/2  •  t/t  o  ] )  and,  after 
to  =  200  ns,  reaches  the  plateau  of  its  constant  maximum 
value.  For  this  current,  the  time  dependent  flow  pressure 
and  temperature  field  and  the  development  of  the  current 
contour  lines  are  numerically  calculated. 

The  basic  equations  are  the  continuity  equation,  the  two 
equations  of  motion  in  the  r  and  z  directions,  and  the  ener¬ 
gy  equation 
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which  are  iteratively  solved  for  a  given  current  contour  by 
an  extended  MacCormack  Code,  as  well  as  the  Maxwell's  equa¬ 
tion  and  Ohm's  Law 


7  X  §  =  no!  (29) 
7  •  B  =  0  (30) 
]  =  a(E  +  (w  x  B)  )  -  -gp-  fj  x  £ )  (31) 

which  yield  the  current  contour  lines  for  a  given  tempera¬ 
ture,  pressure  and  flow  field  by  means  of  a  Relaxation 
Gauss-Seidel  Code  similar  to  that  described  in  Section  2.1. 

The  calculations  proceed  now  in  small  time  steps,  in  which 
for  a  given  current  contour  the  new  pressure,  temperature 
and  flow  fields  are  calculated  by  means  of  an  extended 
MacCormack  Code.  With  these  new  data,  a  new  current  contour 
is  determined  for  the  next  time  step  with  the  new  total  cur¬ 
rent  by  means  of  a  Relaxation  Gauss-Seidel  Code,  and  so  on. 
This  iterative  program  is  quite  extensive  and  needs  the  re¬ 
latively  large  computing  time  of  30-50  min  on  the  Cray  2  at 
the  Computer  Center  of  the  University  of  Stuttgart.  The  re¬ 
sults  so  far  are  quite  promising  and  are  based  on  the  fol¬ 
lowing  initial  values 


3.  Stability  Theor 


3. 1  General  Theory 

The  stability  behavior  of  a  current-carrying  plasma  channel 
has  been  discussed  based  on  a  unique  theoretical  approach. 
One  considers  any  arbitrarily  configuratin  of  such  a  channel 
within  a  fixed  laboratory  system  (xL»yL»zL'  see 
Fig.  23).  The  theory  starts  from  the  equation  of  motion,  the 
continuity  equation  and  the  Maxwell's  equation 

yjr  '  ovj  +  V  { pv^v^  +  p}  =  J*  X  £  (32) 
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After  integration  of  the  equation  of  motion  over  a  channel 
segment  of  length  As  (see  Fig.  23)  one  eventually  obtains 
the  force  balance  for  the  channel  segment  in  the  form 


(  Am) 


X 


AS 


(  38  ) 


where  Am  =  jo  dV  is  the  mass  of  the  plasma 

within  the  segment  of 
volume  AV  and  length  As 

X  As  is  the  sum  of  all  forces  which  act  on  the  segment  and 
accounts  for 

a)  all  gas  dynamic  surface  forces  like  propelling,  friction, 
and  pressure  effects 

b)  the  volume  forces  due  to  electromagnetic  effects  between 
the  current-carrying  elements  inside  and  outside  the  seg¬ 
ment  and 

c)  all  fictitious  forces  like  Coriolis  and  centrifugal  for¬ 
ces  due  to  an  axial  flow  within  a  curved  and/or  twisted 
current-carrying  plasma  channel. 

This  force  X  per  unit  channel  length  has  now  been  calculated 

among  others  as  a  function  of  the  local  channel  geometry 

like  channel  radius  ro,  axis  curvature  R,  channel  twist, 

convergence  or  divergence,  and  also  as  a  function  of  the 

channel  orientation  with  respect  to  an  applied  magnetic 

field,  §o  and/or  flow  field.2  The  transverse  and  axial  com- 
-+■ 

ponent  of  X  follows  by 

2  H.  O.  Schrade,  "Magnetoplasmadynamic  Effects  in  Electric 
Arcs,"  Interim  Scientific  Report,  AFOSR  82-0298,  Aug.  1983. 
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where 

Az  =  Ttro2  channel  cross  section 


ds-  =  ^Ttr°  d¥~  channel  divergence  or  convergence 
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the  first  moments  of  the  current  density  distribution  across 
A z  with  respect  to  the  x-  and  y-axis,  respectively 


I  =  I  u  current  vector  with  u  =  unit  vector  in  the 

2  z-direction  z 

B  the  self  magnetic  induction  due  to  the  current  through 

the  channel  segment 

Bq  magnetic  induction  due  to  the  applied  field 


In  this  calculation  one  assumes  that,  inside  the  discharge 
channel  and  in  its  immediate  vicinity,  the  time  derivation 
of  the  channelbound  impulse  density  can  be  neglected,  or 
that  flow  and  pressure  field  adjust  in  a  relaxation-free 
manner  to  the  3  x  §  forces  and  vice  versa.  This  assumption, 
however,  is  valid  as  long  as  the  change  of  the  channel  dia¬ 
meter  d(2ro)/dt  or  the  motion  of  the  discharge  channel  ) vg | 
is  small  compared  to  the  average  speed  of  sound  within  the 
considered  channel  segment,  i.e. 


d ( 2r  0 ) 


R  - 
-  T 


The  validity  of  this  inequality  can  be  proven  from  case  to 
case  and,  as  discussed  elsewhere,2  is  fulfilled  for  even  ex¬ 
treme  conditions.  Since  one  distinguishes  here  only  between 
stable  and  unstable  channel  configurations  and  does  not  con¬ 
sider  the  real  change  in  time  of  an  unstable  channel,  the 
validity  of  the  above  assumption  or  the  inequality  (^2)  can¬ 
not  be  questioned  any  longer. 

Based  on  the  knowledge  of  X,  one  may  now  determine  "bal¬ 
anced"  and  "stable"  discharge  channel  configuration.  A  bal¬ 
anced  channel  is  characterized  by  the  fact  that  the  trans¬ 
verse  force  zero  along  the  channel  axis.  In  addition, 
such  a  balanced  configuration  is  also  a  stable  one  if  for 
all  possible  small  disturbances  S  off  this  configuration  a 
force  acts  on  the  channel  in  the  opposite  direction  with 
respect  to  £  and  therefore  tends  to  restore  the  original 
balanced  state  (see  Fig.  24). 


Fig.  24  Illustration  of  a  disturbed  axis  configuration. 

The  requirement  for  a  balanced  axis  configuration  is  there¬ 
fore 

Xj_  =  0  (43) 

and  for  a  balanced  and  stable  configuration 

Xj_  •  S  <  0  (44) 

In  the  following  sections  this  criteria  will  be  applied  a) 
to  a  nozzle-type  MPD  arc  in  order  to  investigate  its  stabil¬ 
ity  behavior  and  b)  to  the  spot  conditions  of  the  arc  cath¬ 
ode  attachment  in  order  to  better  understand  the  electrode 
phenomena . 


3.2  Onset  Phenomenon  Explained  b 


a  Flow  Discharge  Instabilit 


N 


expands  out  through  an  opening  or  a  nozzle  throat.  Now,  how¬ 
ever,  instead  of  solid  walls,  we  are  dealing  with  magnetic 
forces  and  an  ohmically  heated  gas  which  should,  however, 
create  the  same  gas-dynamic  propulsion  effect  transverse  to 
the  current-carrying  plasma  tube.  A  typical  banana  or  horse¬ 
shoe-like  plot  of  lines  of  constant  radial  [j  x  £]  forces 
across  a  cross  sectional  area  of  the  bent  tube  is  shown  in 
Fig.  27  and  indicates  a  similarity  to  the  wall  configuration 
of  a  rocket  engine. 
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Fig.  27  Lines  of  constant  radial  [3  x  5]  forces  with 
transverse  plasma  jet. 


According  to  this  electromagnetically  induced  expansion 
flow,  a  propulsion  effect  is  impressed  on  such  a  current 
tube,  or  more  precisely,  a  resultant  gas-dynamic  force  acts 
on  the  control  volume  of  a  tube  segment  (see  Fig.  28).  It  is 

S  =  /  AF  As  +  K ( s+  As )  +  K ( s )  (45) 


where  the  first  term  on  the  right  represents  the  sum  over 
the  surface  forces  around  the  mantel  and  the  two  others  are 
the  forces  which  act  on  both  cross  sectional  endplanes. 


Fig.  28  Illustration  of  the  gas-dynamic  surface  forces 
which  act  on  a  tube  segment  of  length  As. 

It  can  be  shown  that  under  certain  conditions  this  resultant 
gas-dynamic  force  becomes  predominant  and  therefore  bends 
the  current  tube  more  and  more.  As  a  consequence  of  this 
bending,  the  rotational  symmetry  of  the  discharge  will  be 
destroyed,  and  the  nozzle  wall  and/or  anode  ring  on  one  side 
will  be  locally  overheated,  which  leads  eventually  to  an 
anode  spot  and  to  serious  damage  to  the  thruster.  Now  the 
sum  of  all  forces  per  unit  axis  length  given  by  the  surface 
forces,  volume  or  field  forces  and  the  fictitious  forces 
like  centrifugal  and  Coriolis  forces  acting  transversely  on 
a  current  tube  is  given  by  Xj^  according  to  eq.  (39).  Inas¬ 
much  as  one  allows  only  a  bent  and  not  a  twisted  distur¬ 
bance,  the  force  reduces  in  a  first  order  approximation 
of  ro/R  (channel  radius  divided  by  axis  curvature)  to  the 
following  expression 
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The  requirement  eq.  (43)  for  a  balanced  channel  configura¬ 
tion  is  therefore  fulfilled  for  a  straight  discharge  tube 
(R  -*•  ®)  as  shown  in  Fig.  25.  Moreover,  since  any  disturbance 
S  (see  Fig.  26)  points  in  the  direction  of  the  stability 
requirement  follows  by 
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stable 

indifferent  (47) 
not  stable 


It  has  been  shown4  that  instability  should  occur  at  first  in 
a  current  tube  somewhat  downstream  of  the  nozzle  throat. 
This  instability,  however,  affects  then  only  a  part  of  the 
entire  discharge.  Now  the  entire  arc  becomes  unstable  if  the 
conditions  within  the  nozzle  throat  are  such  that  the  apper¬ 
taining  quantities,  as  shown  on  the  left  side  of  eq.  (47), 
become  larger  than  zero. 


One  can  now  model  the  current  density  within  the  nozzle 
throat  in  different  ways  and  obtain  for  the  integral  with 
the  square  of  the  self-magnetic  induction  field  the  ex¬ 

pression 
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where  f  is  a  number  of  the  order  )  ,  which  depends  on  the 
current  density  distribution.  By  taking  a  paraboloid  of 
grade  n,  it  is 


f  =  f ( n)  = 


<n  +  2)2  -  16  £§  +  ^ 


H.  0.  Schrade,  M.  Auweter-Kurt z ,  and  H.  L.  Kurtz,  "Sta¬ 
bility  Problems  in  Magneto  Pi asmadynam ic  Arc  Thrusters," 
AIAA-85- 1633,  18th  Fluid  Dynamics  and  PI asmadynamics  and 
Lasers  Conference,  July  1985. 


For  a  Gaussian  distribution  wit.i  a  cutoff  radius  of  ro  and  a 
half  width  of  rc,  one  obtains 
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The  axial  impulse  flux  in  eq.  (47)  can  now  be  related  to  the 
total  mass  flow  rate  and  the  average  speed  of  sound  within 
the  nozzle  throat.  It  is 
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According  to  Appendix  II  eq.  (11.24),  the  average  tempera¬ 


ture 


Tfc  is  now  a  function  of  the  current  and  the  nozzle 


throat  radius  rt»*  here  one  assumes  a  current  density  dis¬ 
tribution  of  the  shape  of  a  paraboloid.  A  similar  relation 
can  be  obtained  if  one  takes  a  Gaussian  distribution.  In 
that  case,  it  is 
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where  i(ro/rc)  is  now  a  function  of  the  quotient  cuttoff 
radius  ro  divided  by  the  half  width  rc. 

One  now  replaces  in  eq.  (47)  the  two  integral  expressions 
according  to  eqs.  (48)  and  (51)  and  sets  for  Tfc  either  the 
relation  (11.24)  or  (52);  after  some  transformation  one  ob¬ 
tains  the  requirement  for  a  stable  or  unstable  arc  in  the 
form 
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where  now  0  is  a  function  either  of  n  or  of  ro/rc,  depen¬ 
ding  on  which  modelling  of  the  current  density  distribution, 
paraboloid,  or  Gaussian  one  takes.  It  is 
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In  criteria  (53),  the  arc  current  I  is  measured  in  [a],  the 
mass  flow  rate  m  is  measured  in  [kg/s],  M  is  the  molecular 
weight  (e.g.  argon:  =  40)  and  r^  is  the  radius  of  the 
nozzle  throat  in  [m]. 


In  Figs.  29  and  30  the  critical  current  is  now  plotted  as  a 
function  of  the  mass  flow  rate  according  to  eq.  (53)  for  the 
Stuttgart  nozzle-type  thruster  (see  Fig.  2)  with  ru  =  0.012 
and  m  =  0.8  g/s  argon  as  propellant  mass  flow  rate,  taking  n 
and  (ro/rc)  as  parameter,  respectively.  The  curves  with 
n  =  1/2  and  ro/rc  =  2  show  excellent  agreement  with  exper¬ 
imental  onset  values;  both  correspond  to  the  value  4>  =  <t>  ( n ) 
=  <t>(ro/rc)  =  0.31.  In  the  1  g/s  and  1.5  g/s  mass  flow  rate 
tests,  the  onset  conditions  were  not  reached  even  at  5000  A 
and  6000  A,  respectively. 


Based  on  this  analysis,  the  onset  phenomena  are  explained  by 
an  MPD  instability  within  the  flow  discharge  region. 


It  is  a  Known  phenomenon  that  in  a  pulsed  MPD  thruster  ,  , 

and  during  the  starting  phase  in  a  continuous  MPD  thrust- 

8  9 

er  ,  the  arc  cathode  attachment  consists  of  many  hot  spots 
of  high  current  density.  These  spots  are  usually  highly  non¬ 
stationary  and  jump  stat ist ical ly  across  the  cathode  sur¬ 
face.10  In  the  presence  of  an  applied  magnetic  field  paral¬ 
lel  to  the  cathode  surface  these  spots  jump  preferably  in 
the  retrograde  direction,  i.e.  opposite  to  the  Lorentzian 
force  direction.1"  If  the  magnetic  field  lines  are  inclined 
with  respect  to  the  cathode  surface,  one  observes  a  drifting 
of  the  arc  foot  points  or  spots  away  from  the  retrograde  di¬ 
rection.  This  drift  motion  was  experimentally  investigated 
by  Robson"2  and  therefore  is  named  after  him.  This  seemingly 

T.  Yoshikawa,  Y.  Kagaya,  K.  Kuriki,  "Thrust  and  Effi¬ 
ciency  of  New  K-III  MPD  Thruster,"  AIAA  82-1887,  16th  IEPC, 
New  Orleans,  1982. 

°  K.  E.  Clark  and  R.  G.  Jahn,  "Magnetoplasmadynamic 
Thruster  Erosion  Studies,  Phase  I,"  Interim  Report,  AF  Con¬ 
tract  R  0461 1-79-C-0039,  April  1983. 

7  K.  Kuriki  et  al.,  "MPD  Arc  Jet  System  Performance  Test," 
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strange  spot  behavior  can  be  explained  by  a  similar  MPD  in¬ 
stability  shown  for  the  arc  jet  in  the  preceeding  section. 
Now  the  vapor  and  plasma  jet  which  proceeds  from  a  super¬ 
heated  cathode  spot  has  microscopic  dimensions  (of  the  order 
of  microns);  the  jet  is  electrically  conductive  and  thus 
contains  a  current-carrying  channel  (see  Fig.  31)  which 
electrically  connects  the  cathode  spot  and  thus  the  cathode 
with  the  plasma  of  the  interelectrode  gap  and  thereby  with 
the  anode.  To  the  current-carrying  channel  of  such  a  spot 
one  can  now  apply  the  same  balance  and  stability  require¬ 
ments  eqs .  (43)  and  (44)  together  with  eq.  (39)  as  done  be¬ 
fore.  First  of  all,  in  order  to  account  for  not  rotationally 
symmetric  current  density  distributions  within  the  discharge 
channel,  one  models  this  distribution  by  any  "eccentric  pa¬ 
raboloid  of  grade  n".  Such  a  distribution  follows  from  a 
concentric  parabola  through  shifting  the  maximum  to  a  point 
of  the  channel  cross  section  with  the  coordinates 
as  shown  in  Fig.  32.  Based  on  this  modelling  one  obtains2 


J .  5s2  dAz  =  8n  l2  f(n)  f1  +  9(n)  IT  COS0M^ 


VS) 


COS0. 


1  7"  r  u  .  r  r cos 9 1  .  1  -M  M i 

ill  ]  rn  sin9  ^  r  r  n+3  rn  Uin0.J 


where  f(n)  according  to  eq.  (49)  and  g(n)  are  plotted  in 
Fig.  33.  The  transverse  force  per  unit  channel  length 

follows  now  in  third  order  approximation  of  ro/R  from  eq. 
(39)  by 


X1  "  fe  1  t1  +  f<n)  +  R~  (f(n)g(n)  (n+3)  +  1  )ct  +  |  (-£-)2] 


1  r°  dAz 

*Hrs5T  -  /  eV  dAz 

A 

z 


+  R  li  J  0  VVv  +  V<=J  dAJ 


z  y  sy'  z  ‘  r2 


♦&£<*•*>  ['  +ir  ■]-£ 


+  R  fj  /  P  V  fv  +  v  >  dA  }  _  I!  a  ixl, 


Z  X  sx'  z  IR  2 


If  one  now  introduces  the  inclination  angle  5  and  the  dis¬ 
turbance  angle  <0  according  to  Fig.  36  for  such  a  current 
channel  segment  emanating  from  a  fixed  cathode  spot  crater, 
one  can  write  the  transverse  force  per  unit  axis  length  by 
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Fig.  34  Illustration  of  a  disturbed  current-carrying  chan 
nel  at  the  cathode  surface  with  inclined  magnetic  field. 
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where  now 
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k  =  k(n)  =  (n+3)f(n)g(n)  +  1 
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x  I  _  average  self  magnetic  induction  field 

8n  *  roB0  applied  magnetic  induction  field 

(61) 

Herein  one  assumes  that  any  disturbance  t  leads  only  to  a 
bent  (see  Fig.  34)  but  not  to  a  twisted  current-carrying 
channel.  According  to  eqs.  (57)  and  (58),  a  straight  channel 
(with  R  -*■  ®)  is  again  a  balanced  one,  since  =  0.  The  di¬ 
rection  of  any  disturbance  t  points  now  in  the  direction; 
therefore,  a  balanced,  straight  channel  is  also  stable  if 
the  braced  expression  of  the  &/R  component  of  xs  larger 

than  zero.  Neglecting  the  higher  order  terms  of  r0/R,  one 
therefore  obtains  as  stability  requirement 


e  -  -2-  cost,  cos<!>  >  0  (62) 

One  may  now  distinguish  between  the  following  three  cases 

1.  No  applied  magnetic  field,  | B 0  J  ~  ^  =  0  . 

2.  §o  parallel  to  the  cathode  surface,  Z  =  0  . 

3.  B°  inclined  with  respect  to  the  cathode  surface,  Z  *  0. 


Case  1 : 


B, 


-1  =  0 
\  u 


In  this  case  the  stability  requirement  reduces  to 

e  >  0  (stable)  (63) 

which  is  identical  with  that  of  the  MPD  arc  in  section  3.2. 
Here  the  stability  criterion  does  not  depend  on  the  direc¬ 
tion  of  the  disturbance  given  by  $  (see  Fig.  34).  Within 
eq.  (59)  one  may  now  relate  the  average  impulse  transport 

/  pv  2  dA 

A  Z  Z 
z 

along  the  channel  axis  with  the  vapor  and  plasma  pressure 

1  3 

pv  of  the  spot  crater  and  the  ambient  pressure  p^'. 
Therefore,  the  stability  requirement  (63)  can  also  be  ex¬ 
pressed  in  the  form 


Pv  "  P- 


no 

>  3-3  8-s  is 


(64) 


where  I  is  the  current  carried  by  the  spot  and  g  is  the 
s  s 

average  current  density  within  the  spot;  the  factor  3.3  is 
valid  for  a  hemispherical  crater  and  changes  to  about  3.0 
for  a  flat  crater. 

The  pressure  difference  between  the  vapor  pressure  in  the 
spot  and  the  ambient  pressure  must  be  therefore  more  than 
about  three  times  larger  than  the  average  magnetic  pressure 
(PM  =  ^o/8u  I  3  )  in  order  to  maintain  a  straight,  stable 
spot  discharge.  If  this  requirement  is  no  longer  fulfilled 
or  if  e  <  0,  the  channel  becomes  unstable.  Any  small  distur¬ 
bance  will  cause  the  channel  to  bend  more  and  more,  and 
eventually  it  will  contact  the  cathode  surface  at  a  nearby 
site  off  the  original  spot  (see  Fig.  35).  The  new  site  is 
heated  up  and  a  new  spot  is  created.  Again,  for  this  new 
spot,  the  stability  requirement  (63)  or  (64)  must  be  ful¬ 
filled  in  order  to  maintain  a  balanced,  stable  spot.  Since 

1 J  H.  0.  Schrade,  M.  Auwete  r-Kurt  z ,  and  H.  L.  Kurtz, 
"Analysis  of  the  Cathode  Spot  of  Metal  Vapor  Arcs,"  IEEE 
Transactions  on  Plasma  Science,  Vol .  PS-11,  No.  3,  p.  103, 
1983. 


in  this  case  the  stability  criterion  does  not  depend  on  the 
direction  of  the  disturbance,  and  since  one  may  assume  that 
all  possible  directions  of  £  are  equally  probable,  the  spots 
should  move  or  jump  statistically  over  the  cathode  surface. 

Case  2;  B0  parallel  to  the  cathode  surface,  £  =  0  . 

In  this  case,  the  stability  criterion  follows  by 

e  -  y  cos<t>  >  0  (65) 


i.e.  the  channel  stability  depends  now  on  the  direction  of 
the  disturbance  (<j>).  Moreover,  the  force  component  in  the 
tx3  direction  of  eq.  (58)  becomes  unequal  to  zero  in  the 
first  order.  It  is 
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(66) 


For  a  >  0,  i.e.  the  current  density  maximum  is  shifted  to¬ 
wards  the  center  of  curvature,  the  channel  is  less  likely  to 
be  stable  for  any  disturbance  towards  the  retrograde  side 
-n/2  <  <t>  <  +n/2  than  for  a  disturbance  towards  the  amperian 
side  rc/2  <  t>  <  3u/2  (or  -u/2  >  <$>  >  -3n/2)  which  can  be  veri¬ 
fied  by  means  of  criterium  (65).  For  any  disturbed  or  curved 
channel,  there  acts  now  a  force  in  the  ?  x  £  direction  [eq. 
(66)]  such  that  the  curved  channel  always  turns  into  the 
retrograde  direction,  i.e.  £  points  finally  in  the  opposite 
fxS0  direction.  The  channel  therefore  becomes  preferably  un¬ 
stable  by  bending  more  and  more  in  the  retrograde  direction. 


For  a  <  0 ,  the  current  density  maximum  is  shifted  away  from 
the  center  of  curvature  maximum  towards  the  convex  side  of 
the  channel.  According  to  criteria  (65)  and  eq.  (66),  the 
channel  should  now  bend  preferably  in  the  amperian  direc¬ 
tion.  One  can,  however,  show  that  by  accounting  for  the 
higher  order  term  (ro/R)ak  in  eq.  (58),  the  channel  this 
time  assumes  a  curved,  stable  configuration.  The  possible 
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consequence  of  this  is  that  the  channel  does  not  come  in 
contact  with  the  cathode  surface  on  the  amperian  side  of  the 
spot.  This  is  not  the  case  for  a  retrograde  disturbance  with 
a  >  0 ,  for  which  the  channel  is  unstable  and  bends  more  and 
more.  Observable  current  channels  should  therefore  be  ampe- 
rically  bent  or  straight  while  the  observation  of  a  retro¬ 
grade  channel  configuration  would  require  high  time  resolu- 
t ions . 


In  conclusion  to  Case  2,  a  spot  prefers  to  jump  or  move 
either  in  the  retrograde  or  the  amperian  direction,  depend¬ 
ing  on  whether  the  current  density  maximum  is  shifted  to¬ 
wards  the  concave  or  the  convex  side  of  the  discharge  chan¬ 
nel.  Because  of  the  unstable  channel  behavior  during  the 
retrograde  motion  and  the  stable,  bent  channel  behavior  dur¬ 
ing  the  amperian  motion,  one  may  conclude  that  the  average 
jump  distance  is  smaller  and  the  number  of  jumps  more  fre¬ 
quent  for  the  retrograde  moving  spots  than  for  the  amperian 
moving  ones.  For  large  magnetic  fields  ( Me  «  1,  which  means 
about  »  1  Tesla  for  an  experimental  value14  of  I/r 0  107 
A/m)  the  mechanisms  for  amperian  and  retrograde  spot  motion 
are  practically  identical  and  therefore  in  both  directions 
equally  probable. 

Case  3:  £0  inclined  with  respect  to  the  cathode  surface, 
l  *  0  . 

According  to  eq.  (62)  the  current  channel  is  preferably  un¬ 
stable  with  respect  to  a  disturbance  in  the  retrograde  di¬ 
rection  and  assumes  a  bent  but  stable  configuration  for  an 
amperian  disturbance.  The  principle  difference  now  is  given 
through  the  force  component  in  the  ?  x  3  direction,  which  in 
the  first  order  of  ro/R  follows  by 


14  J.  P.  Datlov,  A.  E.  Guile,  and  B.  Jiittner,  "Unipolar 
Arc  Tracks  on  Stainless  Steel,"  Beitrage  aus  Plasmaphysik, 


Vol.  21,  p.  135,  1980. 
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and  which  turns  the  disturbed  or  curved  channel  around  its 
original  straight  axis  either  in  the  retrograde  (a  >  0)  or 
amperian  direction  (a  <  0).  Now,  however,  the  disturbance  o 
finally  does  not  point  exactly  in  these  directions  but  devi¬ 
ates  or  drifts  away  from  them.  The  drift  angle  follows  from 
the  requirement  that  for  this  position  the  turning  force  in 
the  !  x  3  direction  must  be  zero.  From  eq.  (67)  follows 
therefore 


drift 


-  •=—  -T —  tan^ 
2a  ds 


For  a  spot  motion  towards  the  retrograde  side  it  is  a  >  0, 
and  since  the  channel  in  the  direction  of  the  current  is 
supposed  to  converge  (dro/ds  <  0),  the  drift  angle  <)>drift 
is  proportional  to  Z,  as  long  as  the  angles  are  small.  The 
retrograde  moving  spots  drift,  therefore,  under  a  small  an¬ 
gle  with  respect  to  the  proper  retrograde  direction  towards 
the  right  (in  the  direction  of  Bo)  if  the  Bo  field  points 
upward  and  drifts  towards  the  left  (opposite  to  B0)  if  the 
field  vector  points  downward  (see  Fig.  36a, b).  This  spot  be¬ 
havior  has  been  experimentally  observed  and  is  known  as 
Robson  drift  motion.15,16 


1  A.  E.  Robson,  "The  Motion  of  an  Arc  in  a  Magnetic 
Field,"  Procedures,  4th  International  Conference  on  Phenome¬ 
na  in  Ionized  Gases  IIB,  346,  Aug.  1959. 

16  W.  Hintze  and  M.  Laux,  "On  the  Motion  of  Unipolar  Arcs 
in  a  Tokamak,"  Beitrage  Plamaphysik,  Vol .  21,  p.  247,  1981. 
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4. 


Electrode  Investigations 


4 . 1  Theory 

Spot  Initiation:  If  a  plasma  contacts  the  metallic  surface 
of  a  cathode,  a  space  charge  layer  with  a  thickness  of  sev¬ 
eral  Debye  length  is  built  up.  The  metallic  surface  becomes 
electro-negat ively  charged,  while  in  the  plasma  layer  the 
positive  charges  predominate.  The  initiation  of  a  breakdown 
across  this  layer  may  now  be  explained  as  follows:  the  elec¬ 
tric  field  at  the  tip  of  a  protrusion  (see  Fig.  37a)  at  any 
dust  particle  on  the  electrode  surface  becomes  high  enough, 
so  that  by  field  or  thermal  field  emission17  the  current 
density  at  this  tip  or  surface  impurity  increases  to  such  an 
extreme  value  that  a  local  overheating  and  an  explosive 
evaporation  of  cathode  material  takes  place.  Such  a  condi¬ 
tion  requires,  for  instance  for  copper  cathodes,  current 
densities  of  the  order  of  10  A/m  and  more  in  order  to 
overpower  conduction  cooling  within  the  solid  metal  by  elec¬ 
tric  heating. 

As  a  consequence  of  these  high  current  densities,  one  must 
expect  a  strong  ionization  of  the  evaporated  cathode  materi¬ 
al  and  hence  also  an  increase  of  plasma  pressure  in  front  of 
or  around  the  initial  breakdown  point.  Since  this  increased 
plasma  pressure  also  means  a  higher  electron  density,  the 
Debye  length  within  the  high  density  plasma  decreases  and 
therefore  also  the  thickness  of  the  plasma  layer  in  front  of 
the  initial  emission  point.  Due  to  this  effect,  the  voltage 
drop  across  the  entire  plasma  layer  causes  a  higher  electric 
field  at  the  area  where  a  higher  density  plasma  is  present 
than  somewhere  else  where  the  density  is  lower  (see  Fig. 
37b).  The  higher  electric  field,  however,  leads  to  a  higher 

17  G.  Ecker,  "Electrode  Components  of  the  Arc  Discharge," 
in  Erg.  exakt.  Naturwiss.,  Vol.  33,  p. 1 ,  1961. 

18  G.  Ecker,  "Theoretical  Aspects  of  the  Vacuum  Arc,"  in 
"Vacuum  Arcs,  Theory  and  Application,"  ed .  by  J.M.  Lafferty, 
John  Wiley  &  Sons,  New  York,  1980. 
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electron  emission  and  therefore  to  a  concentration  of  the 
electric  current  at  the  initial  breakdown  point.  The  high 
density  plasma  expands  now  into  the  half  space,  thereby 
forming  a  plasma  jet;  simultaneously  f he  high  density  plasma 
is  pressed  onto  the  cathode  surface,  which  melts  and  thereby 
forms  a  crater-like  pit.  Finally,  a  spot-like  discharge  as 
shown  in  Fig.  37c  and  discussed  in  section  3.3  is  initiated. 


Spot  Conditions:  The  high  density,  quasi-neutral  vapor 

plasma  is  now  electrically  connected  with  the  interelectrode 
low  density  plasma  and  interconnected  with  the  overheated, 
molten  inner  crater  surface  by  a  space  charge  layer  with  a 
thickness  of  several  Debye  length.  It  can  be  shown 
the  Debye  length  given  by 
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is  much  smaller  than  the  crater  diameter;  hence  the  plasma 
layer  may  be  treated  by  a  one-dimensional  approach.  Assuming 
different  temperatures  for  the  electrons  Te  and  for  the 
heavy  particles  T  j_  but  taking  Te  and  Tj  constant 
across  the  layer,  the  electric  potential  U  as  a  function  of 
the  distance  x  from  the  molten  crater  surface  has  been  cal¬ 
culated20  and  is  plotted  in  a  dimensionless  form  in  Fig. 
38.  As  parameter,  one  takes  here  the  dimensionless  voltage 
drop,  -(e/kTe)Uo,  between  the  electrode  or  inner  crater 
surface  within  the  range  of  -0.5  and  -10  and  the  quasi¬ 
neutral  plasma  (at  zero  potential). 


H.  O.  Schrade,  M.  Auweter-Kurt  z ,  and  H.  L.  Kurtz, 
"Cathode  Erosion  Studies  on  MPD  Thrusters,"  18th  Interna¬ 
tional  Electric  Propulsion  Conference,  AIAA  85-2019,  Alexan¬ 
dria,  VA,  1985. 

20  H.  0.  Schrade,  "Basic  Processes  of  Plasma  Propulsion," 
Interim  Scientific  Progress  Report,  Grant  AFOSR  82-0298, 
Sept .  1985. 
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Fig.  38  shows  that  the  plasma  layer  width  depends  on 

the  voltage  drop  (-Uo)  and  that  it  can  be  approximated  by 
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where  Mi  is  the  atomic  weight  of  the  evaporated  and  ion¬ 
ized  cathode  material  and  pe  is  the  partial  pressure  of 
the  electrons  in  the  quasi-neutral  plasma  above  the  molten 
crater  surface  given  in  N/m2.  By  introducing  the  ionization 
degree  in  the  extended  form 
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where  no  is  the  number  density  of  the  neutral  particles  and 
ni  is  that  of  the  zi~times  charged  ion  component,  one 
can  replace  the  partial  pressure  of  the  electrons  by  the 
plasma  pressure  p  according  to  the  relation 
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Hence,  the  plasma  layer  width  becomes 
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which  is  for  the  estimated  conditions  of  the  spot  plasma 
( i  £  0.1,  £  50000  K  >  T, ,  U0  <  20  V  and  p  >  0.1  bar  = 


104  N/m2 


T  £  50000  K  >  T . , 

e  8  1 
1  2  •  10~  m,  i.e.  much  smaller  than  the  spot 


size  of  about  1  nm  and  more  in  diameter. 


The  electric  field  Eo  at  the  electrode  surface  becomes  prac¬ 
tically  independent  of  the  voltage  drop  and  therefore  only  a 
function  of  the  plasma  properties  if  eUo/kTe  exceeds  about 
2  (see  Fig.  38).  For  this  quite  reasonable  condition,  the 
field  at  the  smooth,  molten  inner  crater  surface  can  be  cal¬ 
culated  by 


kT 


E(x=0)  =  E0  = 


i  /  m  .  T 
e  1  [/  L  e 
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m  T  . 
e  l 


-  1] 


1/2 
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/  T 

=  2.03  •  10/  [ /m i  ^  -  0.0234] 
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For  a  tungsten  cathode  with  Mj  =  184,  an  ionization  degree 
of  i  =  1  and  a  plasma  pressure  of  p  >  1  bar,  which  is 

realized  in  a  spot  discharge,  one  obtains  for  a  smooth  inner 
crater  surface 

9 


E0, tungsten  ^  ^ 


10 


V/m 
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Since  due  to  these  relatively  large  fields,  small  Taylor 
cones  may  develop  in  the  crater  melt,  and  the  electric 
field  can  locally  exceed  this  value  by  more  than  one  order 

19  2  2 

of  magnitude.  Now  according  to  Good  and  Muller,  the  cur¬ 
rent  density  due  to  thermal  field  emission  is  plotted  in 
Fig.  39  as  a  function  of  the  electric  field  Eo  for  an  elec¬ 
tron  work  function  of  4.5  eV  (tungsten).  It  shows  that 


G.  N.  Fursey,  "Local  Elementary  Processes  at  the  Forma¬ 
tion  and  Operation  of  Vacuum  Discharges,"  presented  at  the 
5th  Workshop  on  Electrode  Phenomena  (Berlin),  May  1982. 

22  R.  H.  Good  Jr.  and  E.  W.  Muller,  "Field  Emission,"  in 
Encyclopedia  of  Physics,  "Electron  Emission  and  Gas  Dis¬ 
charges  I,"  ed .  by  S.  Flugge,  Vol .  XXI,  Springer  Verlag, 
Berlin,  1956. 
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thermal  field  emission  can  explain  the  extremely  high  cur- 

•  2  3 

rent  density  values  which  are  observed  in  spot  discharges 
(up  to  about  10 13  A/m2). 

The  idea  is  that  at  the  tip  of  the  Taylor  cones  the  current 
density  may  assume  values  of  the  order  quite  beyond  that  of 
averaged  microspot  current  densities  by  measuring  the  cra¬ 
ter  sizes  and  appertaining  currents.15  Based  on  measured 
values  of  Is/ro,  i.e.  spot  current  divided  by  spot  radius, 
one  can  now  calculate,  by  means  of  criterion  eq.  (64),  the 
plasma  pressure  which  must  be  maintained  within  the  crater 
in  order  to  keep  the  spot  discharge  stable.  In  Table  1,  the 
difference  between  the  critical  vapor  and  plasma  pressure 
pv  and  the  ambient  pressure  p*,  is  listed  according  to 
the  measured  (ls/ro)  values  for  different  cathode  materi¬ 
als  for  a  flat  and  a  hemispherical  spot  crater  shape. 

Table  1  List  of  (Is/ro)  values  and  the  appertaining  cri¬ 
tical  pressure  differences  (pv  -  p»)  to  maintain  a  sta¬ 

ble  spot  discharge. 


~  [A/m]  (measured) 

(PV  - 

flat 

p  )  .  .  [bar] 

hemispherical 

w 

2.0  •  10  ' 

191 

210 

Cu 

1.4  •  10' 

94 

103 

Mo 

8.3  •  106 

33 

36 

Fe 

( steel ) 

6.0  •  105 

0.  17 

0.  19 

23  S.  Anders,  B.  Jiittner,  H.  Pursch,  P.  Siemroth,  "Inves¬ 
tigations  of  the  Current  Density  in  the  Cathode  Spot  of  a 
Vacuum  Arc,"  Contributed  in  Plasma  Physics,  Vol.  25,  p.  467, 
1985. 


Spot  Dynamics:  The  microspots  are  highly  nonstationary. 
They  propagate  and  disappear  and,  depending  on  the  electrode 
material  (oxide  surface  layer,  surface  chemistry)  surface 
roughness,  overall  temperature  and  type  and  pressure  of  the 
adjacent  residual  gas  or  plasma,  these  microspots  can  clus¬ 
ter  together  to  larger  spots  of  different  sizes.  In  a  trans¬ 
verse  magnetic  field,  these  spots  may  move  d iscont inuously 
or  jump  in  the  retrograde  direction,  i.e.  opposite  to  the 
Lorentzian  rule.  All  these  effects  are  by  no  means  fully  un¬ 
derstood,  but  nevertheless  the  calculations  and  explanations 
in  the  previous  section  3.3  show  that  this  strange  spot  be¬ 
havior  and  even  Robson  drift  motion  can  be  explained  by  a 
magnetoplasmadynam ic  instability  of  the  current  channel 
which  emanates  from  the  cathode  spot.  Based  on  this  analy¬ 
sis,  a  spot  discharge  becomes  unstable  (see  eq.  64)  if  the 
current  Ig  carried  by  a  spot  exceeds  a  critical  value 
given  by 


s,crit 


=  0.55  •  / 


8-rt  ‘ 
HO 


r  o ' 


(PV  -  Pc 
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The  factor  0.55  is  valid  for  a  hemispherical  spot  crater  and 
changes  to  0.58  for  a  flat  one.  Instability  means  that  any 
small  disturbance  leading  to  a  small  bending  of  the  dis¬ 
charge  channel  causes  the  channel  to  bend  more  and  more  un¬ 
til  it  contacts  the  cathode  surface  at  a  nearby  site.  There¬ 
by  the  cathode  is  heated  up,  the  plasma  and  vapor  pressure 
increase  and  since  the  Debye  length  on  the  touchdown  area 
becomes  very  small,  a  new  breakdown  and  second  spot  will  be 
initiated.  Each  of  both  spots  now  carries  a  part  of  the 
electric  current  so  that  both  are  stable.  By  further  in¬ 
creasing  the  current,  at  least  one  of  the  spots  will  split 
again.  The  number  of  spots  will  now  rise  until  the  final 
current  set  by  the  power  supply  is  reached.  The  cathode  sur¬ 
face  will  therefore  be  covered  more  and  more  by  these  micro¬ 
spots,  which  now,  probably  due  to  local  accumulation  of  mol¬ 
ten  electrode  material,  cluster  together  and  build  larger 
spots.  This  spot  motion  should  be  statistical  if  no  external 


magnetic  field  and/or  flow  parallel  to  the  cathode  surface 
is  applied,  since  the  instability  criterion  is  independent 
of  the  direction  of  the  disturbance,  and  since  the  latter  is 
taken  ad  hoc  in  all  directions  equally  probable. 

If  a  transverse  magnetic  field  is  applied,  any  disturbance 
leads  either  to  an  unstable  bending  of  the  discharge  chan¬ 
nel  in  the  retrograde  direction  or  to  a  stable,  bent  channel 
configuration  in  the  amperian  direction.  The  consequence  is 
a  preferred  retrograde  spot  propagation. 

Spot  erosion:  The  erosion  rate  or  material  loss  due  to  a 

spot  discharge  stems  from  two  mechanisms: 
a)  evaporation  of  cathode  meterial  and 
o)  ejection  of  droplets  and  molten  material. 

While  (a)  evaporation  starts  from  the  beginning  of  the  spot 
discharge,  the  second  mechanism  (b)  occurs  after  a  molten 
layer  has  formed  and,  affected  by  the  high  plasma  pressure, 
is  thrown  out  at  the  rim  of  the  spot  crater  and/or  after 
that  if  several  unstable  spots  cause  splashing  on  a  pre¬ 
heated  molten  area.  These  last  mechanisms  need  a  certain 
time  to  develop,  but  they  are  most  destructive  and  cause  not 
only  a  dramatic  material  loss  but  also  a  quite  severe  damage 
to  the  cathode  surface. 
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4 . 2  Experiment 

In  order  to  assess  the  basic  processes  occuring  on  cathodes 
of  MPD  thrusters  and  related  arc  devices,  and  in  order  to 
supplement  the  theoretical  work  on  cathode  phenomena,  the 
following  experimental  arrangement  has  been  set  up  and  tests 
carried  out. 


Experimental  Setup:  The  experimental  arrangement  is  shown 
in  Fig.  40.  Within  a  vacuum  test  vessel,  an  electrode  system 
as  shown  in  Fig.  41  is  charged  with  an  approximately  rectan¬ 
gular  current  pulse  of  a  PFN . 24  It  allows  the  initiation  of 
an  arc  discharge  of  up  to  1  kA  for  about  2  ms  at  pressures 
between  10“  and  10  Pa.  The  anode  consists  of  a  molybdenum 
rod  which  is  rectangularly  oriented  with  respect  to  a  cylin¬ 
drical  cathode.  The  cathode,  a  metal  rod,  is  clamped  between 
two  water-cooled  terminals  such  that  an  auxiliary  electric 
current  can  flow  through  the  cathode  rods  and  so  that  the 
samples  can  be  easily  exchanged.  Using  this  arrangement,  the 
cathode  attachment  can  be  investigated  within  a  defined  mag¬ 
netic  field  and  under  different  temperatures,  depending  on 
the  auxiliary  current  Ia,  on  the  material,  and  on  the  dia¬ 
meter  ds  of  the  samples.  The  transverse  magnetic  induction 
field  at  the  surface  of  the  cathode  amounts  to 


Bu 


^0  la 

n  d 


s 


(76) 


and  the  temperature  of  the  sample  which  now  variates  along 
its  axis  can  again  be  controlled  by  the  auxiliary  current 
setting  and  the  heat  conduction  at  both  clamps.  The  short¬ 
coming  of  this  arrangement,  that  the  magnetic  field  effects 
cannot  be  investigated  independently  of  the  surface  tempera¬ 
ture,  are  overcome  to  a  certain  extent  by  accounting  for  the 
relatively  large  time  to  ohmically  heat  up  or  cool  down  the 

2  4 

P.  Lehniger,  "Designing  and  Building  of  a  Test  Rig  in 
Order  to  Measure  Cathode  Erosion  Caused  by  a  Pulsed  Arc," 
Senior  Thesis,  IRS,  University  of  Stuttgart,  Sept.  1986. 
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samples.  In  order  to  run  the  sample  at  fairly  low  tempera¬ 
tures  with  high  Bo  fields,  one  switches  on  the  auxiliary 
current  simultaneously  with  the  main  current  impulse,  and  in 
order  to  maintain  high  temperatures  with  low  Bq  fields,  one 
switches  it  off  before  starting  the  main  arc  current  pulse. 

The  main  arc  current  pulse  is  initiated  by  means  of  an  igni¬ 
tion  plug  which  ejects  a  plasma  into  the  gap  region  between 
cathode  and  anode  and  hereby  discharges  the  PFN.  A  typical 
cur rent-vol tage  characteristic  of  the  discharge  is  plotted 
in  Fig.  42,  and  the  average  current  I  as  a  function  of  the 
changing  voltage  Ul  is  plotted  in  Fig.  43. 

Measurements  of  the  Arc  Characteristic:  The  measurements  of 
the  U(t)  -  I(t)  characteristic  of  the  main  current  pulse  oc¬ 
cur  according  to  the  following  experimental  set  up  (see  Fig. 
44).  The  arc  current  is  picked  up  by  means  of  a  current 
transformer  (4)  and  its  signal  fed  to  a  transient  recorder 
(5)  (Biomation  Type  802).  The  voltage  pulse  is  measured  via 
an  isolation  amplifier  (7)  and  recorded  by  another  Biomation 
Type  805  (8).  The  pulse  shapes  of  both  the  current  and  volt¬ 
age  traces  are  then  recorded  by  an  x-y  plotter  (9)  and  by  an 
oscillograph  (6)  for  rechecking.  The  pressure  in  the  chamber 
is  taken  by  means  of  a  Penning  vacuum  gauge  in  the  range  be¬ 
low  10' J  mbar  and  above  this  range  by  a  Pirani  pressure 
gauge . 

Calibration  of  the  Cathode  Temperature:  In  order  to  cali¬ 
brate  the  temperature  of  the  cathode  sample  (thoriated  tung¬ 
sten)  as  a  function  of  the  B  or  heating  current  or  as  a 
function  of  the  shunt  voltage  Ush,  the  sample  temperature 
was  measured  at  three  different  locations  with  NiCr-Ni  ther¬ 
mo  elements  up  to  a  temperature  of  almost  1800°  K  and  a  py¬ 
rometer  from  about  1250°  K  to  almost  2500°  K  in  the  middle 
of  the  sample.  The  water-cooling  setting  of  the  mounting 
supports  on  both  ends  of  the  sample  remained  constant,  and 
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the  pressure  within  the  chamber  was  10~  mbar.  The  scheme  of 
the  test  assembly  is  shown  in  Fig.  45.  The  thermocouples 
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Typical  I(t)-  and  U(t)  Characteristics 
during  the  Discharge  Pulse 
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Fig.  44  Test  assembly  to  measure  the  U(t)  and  I(t) 
characteristics  of  the  arc  pulse 
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Cathode  erosion  experiment 
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Fig. 45  Sheme  of  Test  Assembly  for  Temperature 
Calibration 

The  Abbrevations  stand  for: 

TpK  Sample  Temperature 

1^  Heating  Current 

Uh  Heating  Voltage 

Rc,  Shunt  (R  =  2.4‘ICT4  Ohm) 

sn 

Ush  Shunt  Voltage 

Ip^r  Pyrometer  Current 

Ufp-,  Thermoelectrical  Voltage 

1  L 

p  Pressure  in  the  Vacuum-Chamber 

T  Thermocouple  Reference  Temperature 


were  placed  in  small  holes,  and  the  pyrometer  was  oriented 
to  these  holes  in  order  to  approach  black  body  conditions. 
In  Fig.  46,  the  temperature  distribution  along  the  cathode 
sample  is  plotted  for  different  auxiliary  or  heating  cur¬ 
rents. 

Cathode  Measurements:  The  photomicrographs  in  Fig.  47  show 
typical  arc  traces  left  at  the  edge  of  the  main  arc  impact 
area  on  a  thoriated  tungsten  electrode  at  room  temperature. 
The  charging  voltage  and  the  average  current  amounted  in 
this  experiment  to  Ul  =>  160  V  and  about  I  »  550  A,  respec¬ 
tively.  The  lower  picture  is  a  magnification  of  the  lower 
left  crater  in  the  upper  picture.  In  Fig.  48,  the  traces 
left  on  the  main  arc  impact  area  of  a  steel  cathode  are 
shown  for  average  currents  of  a)  a  180  A  and  b)  =»  550  A.  At 
the  lower  current,  one  still  observes  individual  craters  of 
different  sizes  which  partially  overlap  and/or  build  larger 
molten  areas  of  spot  clusters.  At  the  higher  current  (b), 
the  entire  frame  indicates  an  overall  originally  molten 
area,  and  it  seems  reasonable  to  assume  that  the  arc  attach¬ 
ment  was  still  spotty  but  caused  detrimental  "splashing”  of 
the  surface  material.  At  the  edge  of  the  main  impact  area 
(see  Figs.  49a  and  b),  the  molten  zones  are  still  confined 
to  small  individual  craters  of  about  a  few  tens  of  a  micron 
to  several  microns.  In  Fig.  49b,  the  middle  crater  of  Fig. 
49a  is  enlarged.  The  micrographs  in  Fig.  50a, b  show  the  arc 
traces  on  aluminum  within  the  main  impact  area  for  180  and 
550  A  and  look  similar  to  Fig.  48b,  i.e.  like  an  overall 
originally  molten  area. 

It  is  obvious  that  the  arc  attachment  on  a  cold  (room  tem¬ 
perature)  metal  cathode  is  spotty  and  highly  non-station- 
ary.  The  spots  jump  around  the  surface  but  prefer  an  area 
which  is  closest  to  the  anode.  This  area,  heavily  impacted 
by  spots,  starts  to  melt  and  thus  seems  to  become  an  even 
more  preferred  attachment  region.  This  preferred  sticking 
leads  to  an  increased  evaporation  of  cathode  material,  to 
splashing  and  to  droplet  ejection.  In  order  to  avoid  this 
extreme  erosion,  one  must  try  to  obtain  many  small  spots 
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Fig.  46 

Temperature  Distribution  ale 
(Tungsten:  2%  ThC^) 
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Photomicrographs  of  cathode  spot  craters  on  a 
thoriated  (2%  Th02)  tungsten  surface 


a)  Magnification  100  :  1 


b)  Magnification  100  :  1 


Photomicrograph  of  a  steel  cathode  surface  at  the 
main  arc  impact  area  for  an  average  current  of 

a)  i  s  180  A 

b)  I  ~  550  A 


which  rapidly  jump  over  a  large  area  but  which  should  not 
cause  an  overall  molten  surface.  This  last  condition  seems 
to  be  realized  with  a  so-called  thermionically  emitting 

tungsten  cathode  at  fairly  high  temperatures  but  below  melt- 

..19 
mg  point. 

Two  photomicrographs  of  spot  craters  left  on  thoriated  tung¬ 
sten  cathodes  heated  up  to  1500°  K  and  2000°  K  are  shown  in 
Figs.  51a  and  b,  respectively.  At  the  higher  temperature,  an 
alignment  of  the  spots  can  be  seen  parallel  to  the  cathode 
axis,  which  indicates  that  due  to  the  Bu  field  of  the  higher 
auxiliary  heating  current  Ia  (  *  200  A),  the  spots  move  now 
either  in  the  amperian  or--what  is  more  likely — in  the 
retrograde  direction. 

The  erosion  rate  e  has  been  determined  by  weighing  the  elec¬ 
trode  mass  loss  and  by  measuring  the  appropriate  electric 
charge  transferred  across  the  cathode  surface.  In  order  to 
obtain  measurable  quantities  in  mass  loss,  one  adds  up  hun¬ 
dreds  and  even  thousands  of  discharge  pulses.  The  results  of 
these  automatically  run  and  registered  current  pulses  are 
shown  in  Table  2  for  different  materials  at  about  50°  C. 

The  erosion  rates  for  thoriated  tungsten  in  these  tests  are 
by  a  factor  of  about  two  higher  than  those  measured  in  a 
pulsed  nozzle-type  thruster.13  This  discrepancy  is  very 
likely  caused  by  the  different  geometries  of  the  anode.  In 
this  experiment  with  a  pinned  anode  rod,  the  arc  strikes  the 
cathode  in  a  more  confined  main  impact  area,  thereby  causing 
a  clustering  of  microspots  and  a  stronger  droplet  ejection 
than  in  a  coaxial  electrode  device,  in  which  the  spots  are 
spread  over  a  larger  area  and  therefore  do  not  cluster  as 
m  uch . 

Erosion  measurements  on  hot,  glowing  thoriated  tungsten 
cathodes  so  far  did  not  result  in  reliable  and  reproducible 
data.  With  an  increasing  number  of  pulses,  the  brightness  of 
the  glowing  cathode  diminished  noticeably,  despite  a  con¬ 
stant  heating  current.  Very  likely  the  radiation  or  emission 
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Thoriated  Tungsten 
(2%  Th02) 

(not  polished) 


Thoriated  Tungsten 
(2%  Th02) 

{ polished ) 


Pure  Tungsten 
(not  polished) 


Copper 

(not  polished) 


Copper 
( polished ) 


Number  of 

s  u  r  t  d 

Current  Pulses 

Erosion  Rate 

300 

35  t^g/C 

5000 

41  ng/C 

100 

29  tig/C 

500 

34  ng/C 

6000 

42  ng/C 

100 

44  ng/C 

500 

50  ng/C 

5000 

54  yig/C 

200 

73  |ig/C 

800 

77  ng/C 

100 

74  ng/C 

800 

75  pig/C 

constant  of  the  cathode  sample  changed  due  to  convection 
and/or  chemical  effects  of  evaporated  insulator  material  on 
the  tungsten  surface.  Even  several  days  of  outgasing  did  not 
show  a  marked  improvement.  In  order  to  avoid  any  problematic 
insulation  material  and  to  improve  the  mechanism  of  both 
terminal  fastening  mechanisms,  the  design  is  being  somewhat 
altered . 
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6 .  Summary 

In  order  to  gain  a  more  profound  understanding  of  the  funda¬ 
mental  processes  occuring  in  coaxial  arc  jets  and  MPD  arc 
thrusters,  and  in  order  to  eventually  end  up  with  reliable 
design  criteria  and/or  predict  thruster  performances  under 
various  conditions,  the  following  research  program  under  the 
three  headings 

I.  Flow  Discharge  Calculations 

II.  Stability  Theory 

III.  Electrode  Investigations 

has  been  conducted. 


Under  heading  1  ,  the  current  density  distribution,  flow, 
pressure,  and  density  fields  are  calculated  for  a  (quasi-) 
steady  nozzle-type  MPD  thruster,  and  a  computer  program  to 
calculate  the  thrust  as  a  function  of  the  thruster  geometry 
and  as  a  function  of  the  arc  current  and  the  mass  flow  rate 
has  been  written.  These  model  calculations  are  based  on  a 
partially  two-dimensional  approach,  assuming  a  fully  singly 
ionized  plasma  with  an  isothermal  electron  and  an  adiabatic 
ion  behavior.  The  influence  of  the  j  x  §  forces  on  the  axial 
acceleration  and  on  the  radial  density  distribution  is  pre¬ 
sented.  The  results  are  compared  with  experimental  find¬ 
ings.  Under  the  same  heading,  a  rigorous  computer  program  is 
being  developed  which  allows  determination  of  the  pressure, 
flow  and  density  (temperature)  fields  within  a  pulsed,  cy¬ 
lindrical  MPD  accelerator.  As  initial  conditions,  one  as¬ 
sumes  a  homogeneous  radial  current  within  a  confined  axial 
electrode  region  and  a  supersonic  homogeneous  axial  flow. 
For  a  given  current  rise,  the  current  contour  lines  and  the 
flow,  pressure,  and  density  fields  within  the  cylindrical 
flow  discharge  region  are  calculated  as  a  function  of  time. 
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Appendix  I 

Together  with  the  continuity  equation 

7 ( PV  >  =  0  (1.1) 

the  equation  of  motion  (2)  can  be  written  as 

P ( v  •  V ) v  +  7 • P  =  f  x  B  (1.2) 

By  taking  a  one-dimensional  approach  along  the  flow  axis  (  z- 
axis)  and  by  neglecting  friction,  one  obtains 

pv  37  v  +  li  p  =  ^  x  §,z  {I<3) 

where  v  is  the  average  velocity  in  the  z-direction  and  p  the 
average  gas  or  plasma  pressure  across  a  cross  sectional  area 
A  of  the  nozzle.  Now  from  the  continuity  equation  it  follows 
that  the  mass  flow  rate  through  all  areas  A  =  A(  z)  is  the 
same ,  i  . e  . 


m  =  ov.\  =  a  v  Afc  (1.4) 

where  the  index  t  stands  for  the  nozzle  throat.  Applying 
eq.  (1,4)  one  may  now  rearrange  eq.  (1,3)  into 

d  -y-  +  ^  =  —  v  •  ( 1  x  B)  •  A  dz  (1.5) 

m 

In  order  to  relate  the  pressure  p  with  the  mass  density  p, 

l  -  is 


Pe  +  Pi  +  p0  =  Pe  +  PiO 


(1.6) 


P, 


n  kT  =  ?  p  —  T 
e  e  Me 


(1.7) 


art  lal  pressure  of  the  electrons,  Te  is  the  elec- 


j*  ir 
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Hr  j  *%,  **»*•—  « 
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•  tk  ^  t.  ■»  1  “  %  • 

v\»V  « \ 


and  t  is  the  degree  of  ionization  defined 
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by  the  quotient  of  the  electron  number  density  and  the  heavy 
particle  density. 


n0  +  ni 


(1.8) 


The  partial  pressure  due  to  all  heavy  particles  like  ions 
and  neutrals  follows  according  to  their  isentropic  behavior 
by 


p0i  "  p0 


+Pi  *  -S  T 


(1.9) 


where  T  is  the  heavy  particle  temperature,  R  is  the  univer¬ 
sal  gas  constant,  M  is  the  molecular  weight  of  the  plasma 
and  k  is  the  quotient  of  the  specific  heat.  The  index  t  re¬ 
fers  again  to  the  partial  pressure  and  the  mass  density  in 
the  nozzle  throat.  According  to  eq.  (1,6)  together  with 
eqs.  (1,7)  and  (1,9),  it  follows  now  that  for  an  isothermal 
electron  and  an  isentropic  heavy  particle  behavior 


P  *  p  f«  |  t  *  SSiLi  (£_)<->} 

Me  p.  v  p  ' 


(1.10) 


dp  =  {a  ~  T  +  k  (— )K  ^  }  dp 

Me  pt  vpfc' 


(1.11  ) 


Herein  one  assumes  that  the  degree  of  ionization  a  remains 
constant  during  the  expansion  process  (frozen  flow  condi¬ 
tion).  Hence  it  is  also  according  to  eq.  (1,9) 


I  dp  =  d(«  |  Te  In  p  +  ^  S  Tt  (f)1'-'}  (1.12) 


m 


Mb 
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Introduced  into  eq.  (10)  yields 


=  T  /  v<3  x  B)  *A  dz 


(1.13) 


>r  since  at  z  =  zfc ,  one  assumes  a  fully  developed  flow,  i.e. 
^  =  0  and  (]  x  &)  =  0  at  z  =  zfc 


from  eqs.  (1,3)  and  (1,4)  follows 


=  speed  of  sound 


=  /(I  +  a)  J  Tt 


(1.14) 


Herein  one  assumes  that  within  the  throat  the  plasma  behaves 
isothermally .  This  assumption  is  plausible,  since  the  aver¬ 
age  temperature  Tt  of  the  throat  cross  section  is  deter¬ 
mined  by  ohmic  heating  and  radial  heat  (and  possible  radia¬ 
tion)  losses  (see  Appendix  A).  Hence,  according  to  eq. 
(1,14),  the  square  of  the  axial  velocity  at  any  cross  sec¬ 
tion  z  downstream  of  the  nozzle  throat  follows  by 


+T  /  v  (]  x  §)  *A  dz 


(1.15) 


Or  by  dividing  through  vfc  ,  one  obtains  an  integral  equation 

for  (— )  in  the  form 
vt 


(^_)2 

vt 


1  +  (^)|  +  TTa  *  f 1  ^^“'1 


A.  v .  . 

r  t  t  %  <— 1 


_ 

-r~  /  (1  X  B 1  *A(  z )  dz  (1.16) 
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The  thrust  according  to  eq.  (6)  given  by 


Th  ■  1  V  t  Aex  p  +  I  /  3  X  i  dV| 

v2 


(1.17) 


can  now  be  rewritten  by  replacing  pex  according  to 
(1,10) .  It  is 


p  .  p  »,  ,!t 

*ex  Mex  M  t  M  t  1  Pt 


(1.18) 


and  introduced  into  (1,17)  yields 


m 


Th  =  m  vex  +  Aex  pex  |  Tt  (a  +  +  \J  3  x  B  dV  | 


+  -* 


Together  with  the  continuity  equation 


(1.19) 


•K 


“  *  Atptvt  ■  A«pexve* 


and  the  relation  eq.  (1,14) 


V  -  "  +  «>  I  Tt 


one  obtains  for  the  thrust 


Th  = 


=  m  vt  { 


+  1  iL.  [a  +  ( 

V.  1+a  v  L  v 

t  ex 


vt  At  'X-1 

v  A 
ex  ex 


)  ]  +  |  /  )  *  B  dv| 

V2 


(1.20) 
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Calculation  of  the  nozzle  throat  condition 


According  to  the  continuity  equation 


m  =  A  p  v  =  const, 


one  obtains  along  a  flow  line 

dv  +  d£  +  dA  ,  0 
v  p  A 


and  for  the  equation  of  motion  [see  also  (1.3)] 


dv  +  d£_  = 

V  2 

PV 


ti  X  «]. 


dz 


pv‘ 


Since  now  within  the  nozzle  throat 

dA  =  0  and  [3  x  £]z  =  0 

one  obtains  form  eqs.  (II. 1)  and  (II. 2) 
vt2  =  (^)  =  speed  of  sound 


(II.  1 ) 


(II. 2) 


(II. 3) 


and  since  one  may  assume  isothermal  behavior  in  the  nozzle 
throat  according  to  eq.  (1.11)  or  (III. 3),  it  is 

V  -  "  *  “>  Tt  "  "  +  “>  I  Tt 

In  order  to  relate  the  average  temperature  across  the  nozzle 
throat  Tt  =  Te  to  the  arc  current  I,  the  supposition  of 
a  "fully  developed  arc  column"  within  the  throat  will  be 
made,  i.e.  all  the  ohmic  heat  input  per  unit  arc  length  is 
dissipated  away  radially  by  heat  conduction  and  possibly  by 
radiation  (see  Fig.  I I . 1 ) .  This  assumption  is  quite  reason¬ 
able,  since  the  axial  temperature  gradient  is  much  smaller 
than  the  radial  one  and  since  any  heat  input  at  the  smallest 
cross  sectional  area  cannot  increase  the  velocity  beyond  the 
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speed  of  sound.  Hence  the  governing  equation  for  the  nozzle 
throat  is 

■  2 

Vq=l_-RL  (II. 5) 

where  q  =  -  X  T  ur 


is  the  heat  flux  vector  due  to  conduction  in  the  r-direc- 
tion,  and  RL  is  the  radiation  loss. 


For  a  nozzle-type  MPD  arc  thruster,  the  chamber  and  throat 
pressure  is  almost  an  order  below  one  bar,  so  that  the  radi¬ 
ation  losses  are  negligible. 


The  energy  balance 


(II. 5) 

1  £L 

r  dr 


reduces  therefore 

'r  x  If  t>  -  -Ir 


to 


(II. 6) 


or  after  integration  from  0  to  r 


2nr 


(II. 7) 


One  now  takes  the  current  density  as  a  paraboloid  of  the 
n'th  degree 


with 


j  -  j  [i  - 

r  0 


n+2 

n 


ti  r  o ' 


(II. 8  ) 

(II. 9) 


r  o 

where  1=2 n  /  jrdr  is  the  total  current, 

0 


and  the  radii  ro  of  the  conducting  channel  can  be  replaced 
by  the  throat  radii  rf 


With  this  modeling  of  the  current  density  distribution  one 

also  obtains  the  electric  conductivity.  Since  within  the 

throat  f  +  _ 

3=0  E,  wlere  the  electric  field  E 

is  constant  over  the  cross  sectional  area,  one  obtains 


-  100  - 


j(r)  ~  a( r)  =  a  [l  -  (£^)n]  (11.10) 

and  introduced  into  eq.  (II. 7)  yields 

2*r  X  T  =  "  211  r  dr 
or  oO  ro 

=  -  nr 2  {1  -  (11.11) 

a  n+2  r0 


Since  now  within  the  throat  the  electron  and  heavy  particle 
temperatures  are  assumed  identical,  and  since  the  conduction 
losses  of  the  electrons  are  by  far  dominant,  one  may  replace 

2  b 

the  conduction  coefficient  X  by  that  of  the  electron  gas 


k3/2/T" 

X  =  2.08  - - 

/m  Q 
e  ee 


1 


1 


+ 

v 


(11.12) 


which  for  a  singly  full-ionized  plasma  (also  valid  for  not 
fully  ionized  plasmas  if  the  degree  of  ionization  is  larger 
than  ~  30%)  reduces  to 


X 

e 


2.08 


ee 


(11.13) 


Herein,  Qee  is  the  collision  cross  section  between  two 
electrons  and  follows  by 

2  k  3t  3 

-  -  *  !  e  ^2  ^  -,o__2  __3  e 

T  l 


'ee  4  ^4neokT 


In  { 1  +  72  •  it 2  ei 


ne  e 


*  8.77  •  10 


-1°  __L  in  { 5 . 7  •  10-3  -^7*}  [m2] 

Te  pe 


if  Te  is  in  Kelvin  and  pe  in  Pascal.  Hence  it  is 


M.  Schindler,  "The  Transportcoef f icients  of  Plasmas  in 
a  Multicomponent  Approach,"  Diplom  Thesis,  Institute  for 
Plasma  Physics,  University  of  Stuttgart,  1964. 
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k  =  1.27  •  10_i 
e 


T 

In {  5 . 7  •  1 0“  3  Yj^} 


U°K  -I 


(11.14) 


Since  the  In  term  in  the  denominator  is  only  a  slowly  vary¬ 
ing  function  and  for  the  throat  conditions  is  about  3,  the 
coefficient  of  the  electron  heat  flux  is  about  proportional 
to  the  electron  temperature  Te^/2  and  can  be  set  for 
this  problem  to 


V.  =  a  T  5/2  =  4.23  •  1 0- 1  1  T  5/2 
e  e  e 


rWatt  -I 

Lrr>  X  J 


(11.15) 


Introducing  this  into  eq.  (11.11)  yields 


2'c  ‘  TeV2  h  Te  *  '"r2  ¥  11  - 


2  rv  -,n 


Te(r) 

I 


d  T  =  -  ii i  J  {1  -  Jj  (£-)n}  r  dr 
e  2  a  a  r=0  0+2  r° 


and  after  integration,  one  obtains 


<'>  ■  V/2  - *2  [’  -  'hIi'2  (?7)n]  (II-,6> 


The  maximum  value  of  the  electric  conductivity  5  follows  ac¬ 


cording  to  kinetic  theory  by 


e2  n 


Z  n  ✓-  /m  kT  Qq„ 
v  3  it  e  e  ev 
v 


(11.17) 


and  for  a  singly  (more  than  30%)  ionized  plasma,  one  obtains 


3.88  •  10' 


T. 3/2  =  r'T  3/2  &j]  (11.18) 


T 


where  y  is  only  a  slowly  varying  function  of  Tg  and  pg  and 
can  be  set  to 

Y  =  1.30  •  10~5  [ — (11.10) 

VmkJ/z 

Replacing  in  eq.  (11.16)  a  by  eq.  (11.18)  and  j  by  eq. 
(11.10),  one  obtains  for  r  =  ro 


or 


1 


(i 


[T^ 

T 


8u2  ro2  a  Y 


(11.20) 


Since  the  temperature  at  the  edge  Te(ro)  can  be  assumed 
small  compared  to  the  centerline  temperature  Te,  it  fol¬ 
lows  in  good  approximation 


!  -  1 
1 _ n  ]  1/5 

a  y 


(J_)2/5 


(11.21  ) 


According  to  eqs.  (11.10)  and  (11.18),  it  is 

Te(r)  =  Te  [1  -  (|-)n]2/3  (11.22) 


and  therefore  the  average  temperature  Tg  across  the  nozzle 
throat  can  be  related  to  the  maximum  temperature  Te  by 


ru 

J1  [1  -  z"'2  l2'3 

T  = 

2n 

/  T  r  dr 
i,  e 

—  rn 

dz 

1  — 
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e 

itr0 

0 

0 

(11.23) 

Hence,  it 

follows  according 

to  eq 

.  (11.21) 

T  = 
e 

Tt  = 

[1—  J_lV5 
L8u  ay  1 

•  ?(n) 

(i_)2/5 

(11.24) 

with  ¥(n) 

=  (1 

+  1  )V5  . 
n  ' 

1 

J  [1 

0 

-  zn/2  ]2/3  dz 

(11.25) 

plotted  in  Fig.  II. 2. 


Appendix  III 

Calculation  of  the  chamber  pressure  po  based  on  the  nozzle 
throat  conditions  assuming  an  isothermal  behavior  for  elec¬ 
trons  and  heavy  particles  from  the  cathode  chamber  to  the 
nozzle  throat. 

Neglecting  j  x  B  effects,  the  one-dimensional  equation  of 
motion  follows  by 


d(3-)  +  f  -  0 


(III.1 ) 


and  with 


p  -  p  "  k  Te +  o  h  T 

■  0  ('  ♦  “>  Te 


(III. 2) 


and  since  Te  =  Tt  *  const, 
one  obtains 


dp  .  ,,  +  «)  —  Tt  dp 


(III. 3) 


and  the  equation  of  motion  becomes 


v  ^  1c 

•=—  +  (1  +  a)  —  T.  lnp  =  const 
z  m0  t 


(III. 4) 


Since  now  the  velocity  v  =  Vo  in  the  chamber  can  be  set 
equal  to  zero,  one  obtains 


(1  +  a)  lnp0  =  (1  +  a)  T„  lnpk  + 


v.v.v. 
-  ■  W 

«.  v 
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